Electron spin dephasing in two-dimensional systems with anisotropic scattering 



A. V. Poshakinskiy and S. A. Tarascnko 
A.F. Ioffe Physical- Technical Institute, Russian Academy of Sciences, 194021 St. Petersburg, Russia 

We develop a microscopic theory of spin relaxation of a two-dimensional electron gas in quantum 
wells with anisotropic electron scattering. Both precessional and collision-dominated regimes of spin 
dynamics are studied. It is shown that, in quantum wells with noncentrosymmetric scatterers, the 
in-plane and out-of-plane spin components are coupled: spin dephasing of carriers initially polarized 
along the quantum well normal leads to the emergence of an in-plane spin component even in the case 
of isotropic spin-orbit splitting. In the collision-dominated regime, the spin-relaxation-rate tensor 
is expressed in terms of the electric conductivity tensor. We also study the effect of an in-plane 
and out-of-plane external magnetic field on spin dephasing and show that the field dependence of 
electron spin can be very intricate. 
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I. INTRODUCTION 



The spin dynamics of charge carriers in semiconduc- 
tor structures has been attracting a great deal of atten- 
tion [l[ . Much effort is focused on experimental and the- 
oretical studying the electron spin dephasing in quantum 
wells (QWs) and obtaining the controllable spin lifetime 
(for a recent review see Refs.Q-H)- It is established that, 
in a wide range of temperature, carrier density and mo- 
bility, the spin lifetime of a two-dimensional electron gas 
is limited by the D'yakonov-Perel' (DP) spin dephasing 
mechanism [f| 0). The mechanism is based on preces- 
sion of individual electron spins in the Rashba and/or 
Dresselhaus effective magnetic field and is highly sen- 
sitive to the QW crystallographic orientation |7Hlfj|] as 
well as details of electron scattering by structural de- 
fects and phonons. Depending on the ratio between the 
period of spin precession in the effective field and the 
momentum relaxation time of carriers, the spin polariza- 
tion monotonically decays or exhibits damping oscilla- 
tions [IT| - [l3 |. So far, the DP mechanism has been theo- 
retically analyzed for central electron scattering, neglect- 
ing possible anisotropy of scattering potential. However, 
such a model does not always describe electron scattering 
in QWs adequately. Transport measurements reveal that 
electron mobility and scattering rate can be anisotropic 
in the QW plane even for (OOl)-grown structures |15r - 
[ItJ • Strong in-plane anisotropy of electric properties has 
been also demonstrated recently for QW structures with 
embedded semidisk-shaped or elongated dots 0, . 

In the present paper, we study the electron spin de- 
phasing in QW structures with anisotropic scattering 
potential and derive equations for the spin-relaxation- 
rate tensor. We show that anisotropic scattering qualita- 
tively modifies the spin dephasing both in the collision- 
dominated and oscillatory regimes. The paper is orga- 
nized as follows. In Sec. HH we present a general formal- 
ism for describing the electron spin dynamics in quan- 
tum wells in the presence of anisotropic elastic scattering. 
Collision-dominated regime of the DP spin dephasing is 
considered in Sec. IIIII We show that the spin-relaxation- 



rate tensor can be expressed in terms of the constants 
of spin-orbit splitting and the electric conductivity ten- 
sor. In (OOl)-grown QWs with anisotropic scatterers, the 
longest lifetime of electron spin along the growth direc- 
tion is achieved in QWs with structure inversion asymme- 
try where the Rashba effective field is nonzero. The oscil- 
latory regime of spin dephasing, which is realized in high- 
mobility QWs, is considered in Sec. IIVI It is shown that 
anisotropic scattering leads to a coupling between the in- 
plane and out-of-plane spin components even in the case 
of isotropic Rashba or Dresselhaus spin-orbit splitting. 
In particular, the spin dephasing of carriers initially po- 
larized along the QW normal leads to the emergence of 
a net in-plane spin component which then also vanishes. 
We also analyze the effect of an external magnetic field 
on spin dephasing both in the collision-dominated and 
oscillatory regimes and show that the field dependence 
of electron spin can be very intricate. The main results 
of the paper are summarized in Sec. [V] 



II. GENERAL EQUATIONS 

The time evolution of the spin distribution function 
Sfe in the wave vector k space is described by the kinetic 
equation @, 03, HI 



— — + s fe x fi fc = g 

ot 



St s fc , 



(1) 



where f2fc is the Larmor frequency corresponding to the 
effective magnetic field, g is the spin generation rate, 
e.g., due to optical excitation with circularly polarized 
light, and St s& is the collision integral. We consider n- 
doped QW structure with a degenerate two-dimensional 
electron gas and assume that optical excitation creates 
spin polarized electrons directly at the Fermi level, i.e., 
g oc 5(ek — £f), where = h 2 k 2 /(2m*) is the elec- 
tron kinetic energy, m* is the effective mass, and Sf is 
the Fermi energy. Such resonant excitation is commonly 
used in experiments to minimize electron gas heating [12| - 
Hil HH HI! • Under these conditions, the spin dephasing is 
determined by the effective magnetic field and details of 
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electron scattering at the Fermi level, and energy relax- 
ation processes are negligible. For elastic spin-conserving 
scattering, the collision integral has the form [23| 



(2) 



where Wkk' is the rate of electron scattering from the 
state k 1 into the state k and it is assumed that the spin- 
orbit splitting Ml k is much smaller than the Fermi en- 
ergy [24j . Below, we take the scattering rate in the form 
Wkk 1 = 2-kTi 2 / {m* L 2 ) Wkk' S(ek — Sk 1 ) with L 2 being the 
normalization area. 

To solve Eq. (JIJ we decompose the distribution func- 
tion Sfe, the frequency £}k, and the scattering rate Wkk' 
into angular harmonics |20| 



s k = X!' 

n 

n=±l 

wkk- = J2 w ^ einv+im<p ' 



(3) 



where ip = arctan(fcj / /fc x ) and ip' = arctan(fcy/fc^,) are the 
polar angles of k and fe', respectively. The dominant con- 
tribution to the effective magnetic field in quantum wells 
is linear in the wave vector 0, H3- Therefore, we assume 
that the frequency £lk contains only terms with n = ±1; 
the coefficients fl±i are related by fli = The 
coefficients w n ,m satisfy the relations w n>m = w* 



— n,—rm 



It', 



(-1) 



and 



w 



ii.O 



W ,n = (tl / 0) 



which follow from reality of the scattering rate, time in- 
version symmetry, and the optical theorem, respectively. 
Substituting the scries ([3]) for s k , , and Wkk' in Eq. ([1]) 
we obtain the system of linear differential equations for 
the angular harmonics s n 



ds n 
~dt 



m=±l 



s n - m xCl m = gS nfi -w Q , s n + ) w n - m s m 



. (4) 

Here, it is assumed that g is independent of the direc- 
tion of k and, therefore, contains only zero angular har- 
monic. By solving Eqs. (|4]) numerically or analytically 
one can find the time dependence of Sq and thereby the 
evolution of the total spin density S = (l/L 2 ) J^k s fe = 
m*/(2Trh 2 )f™s de. 



III. COLLISION-DOMINATED REGIME 

In this section, we consider the case of frequent electron 
collisions, when the spin rotation angle between scatter- 
ing events is small. In this regime, the anisotropic part 
of the spin distribution function s k is much smaller than 
Sq and Eqs. ((H) can be solved iteratively [1, 0]- Such a 



procedure gives the following equation for the spin den- 
sity 



f = G-TS. 

dt 



(5) 



where G = (1/L 2 ) J2k 9 i s the total spin generation rate 
per unit area and T is the spin-relaxation-rate tensor. 
The latter is defined by 

TS= {l/L 2 )Y^{s-i x Oi +si x n_i) , (6) 



where, to first odcr in the effective magnetic field, the 
harmonics a±i are to be found from the equation 



s x ft k — St s k 



(7) 



The calculation of Eqs. © and ([7]) is similar to the 
calculation of an electric current density j induced by a 
static electric field E. Indeed, the current density is ex- 
pressed via the electron distribution function f k by j = 
(e/L 2 ) ^2 k (vif-i+v-ifi), where e is the electron charge, 
v±i and f±i are the angular harmonics of the electron 
velocity v k = hk/m* and the distribution function, re- 
spectively Within linear in E regime, the harmonics f±i 
are found from the equation e(dfo/dsk)v k ■ E = St/fc, 
which is similar to Eq. ([7]). Such an analogy allows us 
to express the spin-relaxation-rate tensor in terms of the 
tensor 2 x 2 of in-plane electric conductivity er as follows 



7TTO 



[l 3 Tr(AerA T ) - A<tA t ] 



(8) 



where J 3 is the unit matrix 3 x 3, A is the matrix 3x2 
relating components of the frequency tt k and the wave 



vector fe, flk = Afc, and we used that er 



Equa- 



tions (O and © describe the spin dynamics of a degen- 
erate two-dimensional electron gas for arbitrary elastic 
scattering and generalize previous results. If the scatter- 
ing potential is central, then the conductivity tensor is 
diagonal and can be expressed via the momentum relax- 
ation time t\ at the Fermi energy by er = (7/2, where 
er = Tie 2 kp j (2-Km*) , kp is the Fermi wave vector, and I2 
is the unit matrix 2x2. In this particular case, Eq. ((5J has 
the form T = (nk%/2) [j 3 Tr(AA T ) - AA T ] , in agree- 
ment with the result of D'yakonov and Kachorovskii 0]. 

To analyze Eq. © in more detail we consider QW 
grown along z || [001] crystallographic direction. In such 
structures, the matrix A has nonzero components 



A XJ/ = a + (3 , A yx = j3 — a, 



(9) 



where a and (3 are the constants of the Rashba and Dres- 
selhaus spin-orbit splitting, respectively, x || [110] and 
y || [110] are the in-plane axes [a [HI- Then, components 
of the tensor T take the form 



xx — o 



(a - p) a xx , T yy = —^-(a + P) a yy , 



xy " r y X — 2 ipt $ ^)&xy j ^ zz — ^ xx H~ ^yy • (^^) 
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Dependences of the tensor components T xx , T yy , and is symmetric, i.e., er T (0) = er(0). In the presence of 
T xy on the ratio a/(3 for QWs with strong scattering magnetic field, the tensor ct{uj c ) contains both symmet- 



anisotropy are plotted in Fig. 1 by dashed curves. The 
in-plane eigen values 71 and 72 of the tensor T and the 
out-of-plane value j z , which coincides with T zz , are found 
from the equation det(7/3 — T) = and presented in 
Fig. [T]by solid curves. One can see that the minimum of 
j z , which corresponds to the longest lifetime of the spin 
component S z , is achieved in asymmetric QWs where the 
Rashba constant a^O. This is in contrast to QW struc- 
tures with central electron scattering where the Rashba 
splitting is known to decrease the spin lifetime. The anal- 
ysis of Eqs. (fT0|) shows that, at fixed j3, the rate j z reaches 



the minimum at a/ (3 = (a xx — &yy)/ Tr er. We also note 
that 71 7^ 72 no matter how small the ratio a/ (3 is if the 
eigen axes of the conductivity tensor do not coincide with 
x and y, i.e., a xy ^ 0. 

Now we consider the effect of an external magnetic 
field B on spin dephasing. The magnetic field causes 
the Larmor precession of electron spins and cyclotron 
motion of electrons in QW plane with the frequencies 
flz = gf-isB/h and uj c = eB z /(m*c), respectively [27], 
l28j |. Here, g is the effective electron (/-factor, [Ib is the 
Bohr magneton, e is the electron charge, and c is the 
speed of light. Both effects are theoretically described in 
the framework of kinetic approach, with kinetic equation 
having the form [Io|, H3, Hi] 



ds k 

at 



S k X (fife + fix) - UJ C 



k x 



d_ 
dk~ 



Sk 



St Sfc 



(11) 

Solution of Eq. (fTT|) shows that, in the collision- 
dominated regime, the time evolution of the spin density 
S is described by Eqs. and © where (i) the additional 
term S x fi^ is added to the left-hand side of Eq. (J5J and 
(ii) er in Eq. ([5J is replaced by the transposed tensor of 
electric conductivity in the magnetic field er T (cj c ). We 
note that in zero magnetic field the conductivity tensor 
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FIG. 1: (Color online) Dependences of the spin-relaxation- 
rate tensor components T xm , F yy , and T xy (dashed curves) 
and the tensor eigen values 71, 72, and 7 Z = T zz (solid curves) 
on the ratio ce//3 calculated for a xx /Trcr — a xy /Trcr = 1/4. 
The curves are normalized by To = (7ml* /e 2 ) f3 2 Tr er. 



ric cr 

«■(«) 



;s) K) 

er(w c 



= [<t(oj c ) + er (w c )]/2 and antisymmetric 
er T (w c )]/2 parts. Accordingly, the right- 
hand side of Eq. (|SJ) can be also reduced to the sum of 
symmetric T^ s > (w c ) and antisymmetric r^ a ' (w c ) tensors. 
The symmetric tensor rw(u; c ) describes spin relaxation. 
The antisymmetric third-rank tensor rW(w c ) is equiv- 
alent to a pseudovector Jflz, and represents, in fact, a 
correction to the Larmor frequency 

IS EES- There- 
fore, the equation describing the time evolution of spin 
density in the magnetic field has the final form 

A C 

— + S x (fix, + 5il L ) = g - T(lu c )S , (12) 



where 

r 



( Wo ) = ^- {/ 3 Tr [A*^(u c )A T ] - Acr^i^} 

(13) 

is the spin-relaxation-rate tensor, 



2e 2 ^ 



"a/37 



Acr( a )(a; c )A 2 



01 



(14) 



is the Larmor frequency correction caused by cyclotron 
motion, and e a( g 7 is the antisymmetric third-rank tensor 
(Levi-Civita symbol). 

As follows from Eq. (|12[) . the precession of total elec- 
tron spin is determined by the frequency Ol + SQl- The 
frequency correction (Sfi^ depends on the magnetic field 
non-monotonically: It is proportional to the magnetic 
field at small fields, reaches maximum at u c t\ ~ 1, and 
decreases with the further field increase. It may hap- 
pen that fix and <5fi_L have opposite signs and com- 
pensate each other at a certain magnetic field. Such a 
compensation results in a peculiarity in the magnetic 
field dependence of the electron spin. As an exam- 
ple, we consider the simple case of continuous spin gen- 
eration, central electron scattering, and the magnetic 
field B pointed along the QW normal [001]. Then, the 
conductivity-tensor components have the form <j xx (lo c ) = 



'mi 



(w c ) = er/[l + (cj c ti) 2 ], a xy (u> c ) 



-<Jy X (U} C ) 



au> c Ti/[l + (cl> c ti) 2 ], and the spin-relaxation-rate tensor 
(fT"3"|) is diagonal in the chosen coordinate frame (x,y,z). 
Straightforward calculation shows that the components 
of the steady-state spin density S have the form 

Tyy(UJ C ) G X - (fiL + SflL)Gy 
T XX (UJ C ) Tyy(u} C ) + (fi£ + <5fii) 2 

T xx (uj c ) Gy + (fiz, + 5Ql)G x 
T xx (u> c ) T yy {io c ) + (fi L + 8Vl l ) 2 ' 

G z 



S 11 



r zz (w 



(15) 



where 



r(w c ) = r(0)/[l + Kn) 2 ], T xx = Tl k 2 F (a~l3) 2 /2 
Tyy = Tl k 2 F {a + f3) 2 /2, T zz = T xx + T yy @, and 



SSl L = 



Tik 2 F {a 2 - f3 2 ) uj c t x 

2 1 + (w c n) 2 



(16) 
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Shown in Fig.[2]are the magnetic field dependences of the 
in-plane components S x and S y calculated for the spin 
generation G ]| ar, the Rashba spin-orbit splitting, and 
ili/uj c = ±0.01. Such ratios of the Larmor to cyclotron 
frequency can be realized, e.g., in GaAs/AlGaAs QW 
structures [3l|. The dependences plotted for VLl/ui c = 
0.01 (dashed curves) are rather simple: S x is maximal at 
B = and monotonically decays with the field increase; dt 
S y oc B at small fields, reaches maximum, and then de- 
cays. In contrast, the magnetic field dependences of S x 
and S v calculated for Sl L /u c = —0.01 (solid curves) are > 



completely different. The component S x reaches maxi- 
mum at a finite magnetic field corresponding to lo c ti « 2. 
S y has two extrema for a fixed direction of B and changes 
the sign approximately at the magnetic field where S x is 
maximal. Such a behavior is caused by interference of 
VLl and 8£Il which compensate each other at uj c t\ k. 2 
for the parameters chosen. The vanishing of the total 
Larmor frequency leads to the Hanle-like curves in the 
vicinity of this magnetic field, see Fig. [2] The fact that 
S x in the point of compensation is much larger than S x (0) 
is caused by a slowdown of the DP spin dephasing by cy- 
clotron motion [27|, see Eq. (TTSl) . 




FIG. 2: (Color online) Magnetic field dependences of the spin 
components S x and S y calculated for G || x, B || z, the 
Rashba spin-orbit splitting, akFTi = 0.3, and SIl/^c = —0.01 
(solid curves) or Q,l/lo c = +0.01 (dashed curves). 



IV. OSCILLATORY REGIME 

Now we consider the oscillatory regime of spin de- 
phasing which occurs if the relaxation time T\ is longer 
than 1/fife pdl4l4l | . For arbitrary fife and scattering rate 
w(<p, f'), Eqs. (j4]) can be solved only numerically. There- 
fore, we focus on scattering anisotropy and assume the 
isotropic spin-orbit splitting of the Rashba type. The 
scattering rate can be conveniently presented as the sum 
of two terms [23| 

w(<p, ip') = w^(\p - tp'\) + 8w{<p, V ') , (17) 

where w^{4>) = f Q * w(ip' + 4>, ip') dip'/ (2tt), Sw(tp, <p') — 
w(ip, ip') — w^ (\ip — f'\), and we assume that Sw <C w^ c \ 



The term w^ c \<j)) describes central scattering. The cor- 
responding collision integral is expressed in terms of the 
relaxation times r n of angular harmonics of the distri- 
bution function, l/r„ = f 27r w^(cj))(l — cosncj)) dc/>/(27r). 
Then, Eqs. (0} take the form 



ds. 



n 



S-,n-l + S+,n+l 



r- °z,n 
g Z n ,0 

(18) 



dt 



s z , nT i — g±S n _o — 



(19) 



where s ±i „ = s X}7l ± is y ^ n , g± = g x ± ig y , and fl R = ak F 
is the precession frequency corresponding to the Rashba 
field at the Fermi level. 

In the regime of continuous spin generation, when 
g is independent of time, Eqs. (|18p and (fT9"| can be 
solved iteratively. To first oder in Sw, spin compo- 
nents depend on the harmonics Swu = 8w*L 1 _ ± and 
(5u>2,-i = 8w*_ 2 1 = —Suu—i^ = —5u'l _ 2 . Such angular 
harmonics in the scattering rate reduce space symmetry 
of the system and, in fact, correspond to a symmetric 
tensor and an in-plane vector, respectively. Accordingly, 
we define 2x2 tensor Q by Q xx = —Q yy = HeSwi^i 
and Q xy = Q yx = — Imdwi^i and the vector u by 
v x = Re<5w2,-i and v y = — Im 8w2,-i- In these nota- 
tions, the in-plane S\\ = (S x ,S y ) and out-of-plane S z 
components of the steady-state spin density are given by 



S z 



T2 



G z 



l R T l 
T 2 V 



G \\ ~ 7yrQ G \\ 
Gii 



r 2 v 
~&R~ 



9., 



(20) 
(21) 



where G\\ = (G x , G y ) is the projection of generation rate 
onto the QW plane. Equations (|2"0"|) and (|2ip show that 
the in-plane and out-of plane spin components are cou- 
pled in structures with anisotropic scattering; the cou- 
pling strength is proportional to v. In particular, the 
generation of electron spin along the QW normal, i.e., 
G || z, leads not only to S z but also to Sn oc vG z . 
Moreover, even in the case of small scattering anisotropy, 
the in-plane and out-of-plane spin components can be 
comparable to each other provided 91rJ\ is large enough. 
The second term on the right-hand side of Eq. (|2U|) de- 
scribes the in-plane anisotropy of spin dephasing due to 
anisotropic conductivity which, to first order in Sw, has 
the form er = Tie 2 k 2 F /(2TTm*)(I 2 + t x Q). The effect of 
conductivity anisotropy on spin relaxation was consid- 
ered in Sec. IIIII Below we focus on the coupling between 
S\\ and S z and assume, for simplicity, that Q = 0. 

The coupling between the in-plane and out-of-plane 
components of the spin density can be also studied in 
experiments with high time resolution. Shown in Fig. [3] 
are the time dependences S z (t) and S x (t) after a short 




FIG. 3: (Color online) Time dependences S z (t) and S x (t) 
after a short optical pulse orienting electron spins along z. 
The curves are calculated for v || x, vt\ = 0.1, r n = n, and 
two different fi ijn. 
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FIG. 4: (Color online) Microscopic mechanism of the genera- 
tion of in-plane spin polarization when electron spins are ini- 
tially oriented along the QW normal. Precession of electron 
spins in the effective magnetic field followed by anisotropic 
electron scattering and subsequent spin precession in the ef- 
fective field results in a spin polarization along x. 




FIG. 5: (Color online) (a) Example of a noncentrosymmet- 
ric scatterer. Disk with different edges which diffusively and 
specularly scatter electrons, (b) Electron trajectories in QW 
plane. Trajectories shown by solid and dashed lines are con- 
nected by space inversion. 



circularly polarized optical pulse which orients electron 
spins along z at t = 0. The curves are obtained by solving 
Eqs. (|20|) and (|21[) numerically for nonccntrosymmctric 
scattering potential with v || x. One can see that S z (t) 
demonstrates damping oscillations as is expected for the 
oscillatory regime of spin dephasing [ll| . The oscillations 
are caused by precession of individual electron spins in 
the effective magnetic field. The in-plane spin component 
S x is zero right after the pulse, emerges at the time scale 
of momentum relaxation, and then also decays. For the 
parameters given in caption to Fig. [3J S x (t) reaches a few 
percent of S z (0). 

Microscopic mechanism of the generation of the in- 
plane component S x is a three-stage process illustrated 
in Fig. |4j At the first stage [Fig.[4ja)], electron spins ini- 
tially oriented along z precess in the Rashba field with the 
frequency fife. The precession forms a spin distribution 
function Sfc containing the first angular harmonic. The 
electron scattering by noncentrosymmetric defects mod- 
ifies Sfc and, due to the terms oc 8w2 t -i and oc 8w _2,i 
in the collision integral, partially transforms the first an- 
gular harmonic into the second harmonic. The spin dis- 
tribution function described by the second angular har- 
monic is shown in Fig.SJb). Finally [Fig.|4jc)], the subse- 
quent rotation of electron spins in the Rashba field results 
in a net spin polarization of carriers along the x axis. 

The coupling between S|| and S z can be also under- 
stood by analyzing electron trajectories in QW struc- 
tures where the angular dependence of scattering rate 
contains the harmonics 5w2,— i and Sw—2,i- An example 
of scatterer providing such harmonics is a disk, one edge 
of which reflects electrons specularly while the other one 



scatters electrons diffusively, see Fig. Eta). Obviously, 
such scatterers with a preferred orientation in the QW 
plane break the in-plane space inversion [33| . The scat- 
tering anisotropy modify electron trajectories, which af- 
fects the spin dynamics. Indeed, in QWs with centrosym- 
metric scatterers, each electron trajectory has on average 
its counterpart connected by space inversion, see solid 
and dashed lines in Fig. [5jb). Electrons the spin Sq ini- 
tially oriented along z move in the QW plane and gain in- 
plane spin components sy due to rotation in the Rashba 
field. However, the particles propagating along the paths 
interconnected by space inversion gain the opposite pro- 
jections S|| leading to a vanishing average in-plane spin 
polarization. In QWs with noncentrosymmetric scatter- 
ers, the space-inversion symmetry of electron trajectories 
is broken, which results in a non-zero in-plane spin po- 
larization. To first order in scattering anisotropy, Sy is 
determined by the angular harmonics 8w2.-i and <5w_2,r- 
Other harmonics Sw n . m describing noncentrosymmetric 
scattering can also couple the in-plane and out-of-plane 
spin components in higher orders in 5w. 

Equations (JTHJ)- (EU) ar< 3 obtained for the Rashba spin- 
orbit splitting from general Eq. Similar calculations 
can be carried out for the case of Dresselhaus splitting. 
One can see that Eq. ^ is invariant to the replacement 
of the Rashba field with the Dresselhaus one, which differ- 
in sign of fly , and the simultaneous inversion of s x and g x 
signs. Thus, Eqs. ([20]) and (f2"Tj) where FIr, S x , and G x are 
replaced by fljj = ftkp 7 —S x , and —G x , respectively, de- 
scribe the steady-state spin density in (OOl)-grown QWs 
with the Dresselhaus splitting. In general case, if both 
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the Rashba and Drcssclhaus contributions to the effective 
field are present, a prerequisite for the coupling between 
S\\ and S z remains the lack of inversion symmetry in scat- 
tering potential. We also note that the coupling is absent 
irrespective of the form of w{ip,ip') if |a| = |/3| because, 
in this particular case, the frequency Q.k depends only 
upon one component of the wave vector. 

Finally, we discuss the effect of an external magnetic 
field on spin dynamics in the oscillatory regime. Equa- 
tions (|18|) and (fT9|) with the Larmor precession and cy- 
clotron motion being taken into account have the form 

ds zn S_|_ n +l + S— ,n—l , ."L,-S+, n — ^L, + S-.n 

Mr n 1- 1 — o 



the external in-plane magnetic field equal to the effec 



g z S n . - (— - inujc) s 2 ,„ + Sw. 



J n,—m J z,m 3 



■ + VIr S 2i „ t i ± iflL,± S z ,n = g±S n ,0 

( inu c T i^L.z ) s±. n + 



n,— m °±,m i 



(22) 



(23) 



(7/ 



dt 



where Ql,± — ^l.x H^L.y Equations (1221) and ([23)1 are 
valid for arbitrary strength of spin-orbit splitting VLRT n 
and angular dependence of the scattering rate. 

Dependences of the steady-state components S z and 
S x on the in-plane magnetic field B \\ y for continuous 
spin generation along the z axis are shown in Fig. [5] The 
curves are obtained by solving Eqs. (f2"2")l and (f2"3")l numer- 
ically for is || x and different SIrT\. One can see that the 
curves drastically depend on the parameter £IrT\. At 
£Irti = 1, the dependences S Z (B) and S X (B) are simi- 
lar to conventional Hanle curves. The only difference is 
that S x (0) 7^ due to scattering anisotropy, see Eq. (f20| . 
The dependences S Z {B) and S X {B) calculated for large 
UrTi, UrTi = 3 in Fig. [HI look completely different. In- 
stead of a monotonic decrease with the magnetic field, 
S Z (B) increases first with the field, reach maximum at 
Ql ~ CIr, and then decreases. S X (B) is nearly indepen- 
dent of B at small magnetic fields and exhibits a sharp 
rise at w Qr. Such a behavior is caused by a par- 
tial suppression of the DP spin dephasing mechanism by 
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FIG. 6: (Color online) Dependences S z and on the in-plane 
magnetic field B \\ y measured in units of Q.l/£Ir- The curves 
are calculated for G \\ z, v || x, r n = ri, vt\ = 0.1, and two 
different Qrti. 



tive field in high- mobility structures [32(. We also note 
that the dependence S Z (B) is always even despite the 
fact that the vector S(0) is not aligned along the z axis. 
The evenness of S Z (B) follows from Eqs. ([22]) and (j23)) . 

For the magnetic field pointed along the QW normal 
(f2^ || z) and continuous spin generation, the steady- 
state solution of Eqs. ([22]) and (|23|) can be found analyt- 
ically. To first order in the scattering asymmetry 5w and 
for G II z, the solution has the form 



s v 



Re 

Re 
1 4 



(i^ - £l R fif 2 



1 + fflif2 + 2iQ L (l/n - iw c )/fi|, 
+ w x ) fif 2 



1 + if2 L f 2 + 2iOi(l/ri - iu c )/VL 2 R 
(uj c + & l ) 2 tI 



-G z 



(24) 



where l/f„ = 1/t„ — mw c — iSl^. The z-component of 
the spin density quadratically increases with the mag- 
netic field growth, which is caused by a slowdown of the 
D'yakonov-Perel' spin dephasing mechanism in the per- 
pendicular magnetic field [H, [29| ■ The dependences of 
S x and S y on the magnetic field B \\ z are more com- 
plicated and drastically depend on parameters. Exam- 
ples of such dependences are shown in Fig. [7] First, we 
note that S x ^ in zero magnetic field due to scattering 
anisotropy. With the field B increase, S x decreases and 
changes the sign. The component S y depends linearly 
on B at small magnetic fields, reaches an cxtremum, and 
then decreases. The curves calculated for Qz,/cj c = —0.05 
(solid curves) have additional Hanle-like peculiarities at 
uj c t\ » 3.2. The analysis of Eqs. (f2"4")l shows that the 
peculiarities occur at lo c « Q.r^J—uj,./ (2Q,l) and are of 
similar origin as those in Fig. [2] caused by nulling the 
total frequency fix + SQl. 



0.1 



. N 



0.0 



-0.1 - 



s 

G 


z 

AS 




/ 


y\ . y j 




/ *S 








>^ A 


\s / \ 




L v 








s \ I 



FIG. 7: (Color online) Dependences S x and S y on the mag- 
netic field B || z measured in units of uj c ti- The curves are 
calculated after Eqs. (|24l) for v || x, r n — n, vr\ = 0.1, 
Q,rti — 1, and SIl/u c — —0.05 (solid curves) or Ql/u c = 
+0.05 (dashed curves). 
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V. SUMMARY 

We have developed the microscopic theory of elec- 
tron spin dcphasing in QW structures with anisotropic 
scatterers. It is shown that, in the collision-dominated 
regime of spin dcphasing, the spin-relaxation-rate tensor 
is determined by constants of spin-orbit splitting and the 
electric conductivity tensor. In (OOl)-grown structures 
with anisotropic in-plane conductivity, the longest spin 
lifetime of electrons polarized along the QW normal is 
achieved in asymmetric QWs with a finite Rashba split- 
ting. We have demonstrated that, in structures with non- 



ccntrosymmctric scattering potentials, the in-plane and 
out-of-plane spin components are coupled to each other. 
The coupling is caused by breaking the space-inversion 
symmetry of electron trajectories in the QW plane and 
is more pronounced in structures with strong spin-orbit 
splitting. The engineering of electron trajectories pro- 
vides an additional approach to manipulate electron spins 
in low-dimensional semiconductors. 
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